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This paper considers the effect of electrostatics on the stability of a charged membrane. We 
show that at low ionic strength and high surface charge density, repulsion between charges on 
the membrane renders it unstable to the formation of holes. A free straight edge is unstable to 
modulations with wavelength longer than the Debye screening length. Hence at low ionic strength, 
membranes will disintegrate into vesicles. We use these results to interpret stable holes in red blood 
cell ghosts (Steck et. ai, Science 168:255(1970)). 
Lipid membranes, the sheets which form the boundaries of cells, must maintain their integrity for cells to live 
0. Since a tear in a membrane exposes the hydrophobic interior of the sheet to water, tradition maintains that 
membranes have a high line tension (energy per unit length of exposed edge) and never spontaneously form holes. We 
show that membranes can form holes in a simple model of a charged membrane. By examining the competition between 
electrostatic repulsion and line tension, we find a parameter range where electrostatics dominates and membranes are 
unstable to hole formation. Although much previous work examines electrostatic effects on lipid assemblies Q and 
the bending moduli of membranes P|3,p[ , we are only aware of one study addressing electrostatics and line tension 



FIG. 1. Image of hole in red blood cell ghost, from Lieber and Steck |g] Figure 7. 

These calculations were motivated by experiments which observed stable holes in red blood cell ghosts |^J^,p|,^0[ , 
the size of which depends on the ionic strength of the surrounding fluid |f| (Fig. 1). We compare our calculation with 
these experiments, and argue that the observed morphologies require the spectrin skeleton. 

Consider a thin, axisymmetric membrane with constant charge density a on each side (Fig. 2a). The membrane 
attracts ions of the opposite sign, creating a screening layer (of thickness k _1 in Debye Huckel theory). We assume 
that the membrane is flat; this is a good approximation when the hole radius R and screening length are small 
compared to the membrane radius of curvature. The plane of the membrane is at z = 0. Throughout this paper, the 
membrane area and charge are held constant. 
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FIG. 2. (a) Sketch of a membrane hole; shading denotes constant charge density, (b) Sketch of the screening cloud (shaded) 
around the membrane (bold line). When R < n~ x (top), the screening cloud fills the hole. When R > K _1 (bottom), the cloud 
extends only k -1 from the edge. 

Repulsion between charges on the membrane favors the creation and expansion of holes. Heuristically, the electro- 
static energy is determined by the volume of the screening cloud: the total charge in the screening cloud is constant, 
so expanding the cloud lowers its energy. First, when the hole is small relative to a screening length (Fig. 2b), the 
cloud expands by a volume ~ R 2 k~ x over the hole. Since the electric field strength is ~ cre^ 1 (e w is the dielectric 
constant of water), the energy decrease is ~ a 2 n~ 1 e~ 1 R 2 . The energy increases because an edge of membrane is 
exposed; if the membrane has a line tension 7, the energy cost of a small hole is U ~ jR — a 2 n~ 1 e~ 1 R 2 . For small R 
line tension dominates, and a hole closes. 

In the opposite limit, the hole is much larger than a screening length. The cloud does not fill the hole, but inhabits 
an additional volume equal to that of a tube of radius k^ 1 and length R. The total energy is U ~ (7 — a 2 n~ 2 e~ 1 )R. 
Depending on 7, a, and k, the energy can be positive or negative. The latter case corresponds to instability and hole 
growth. 

We quantify this argument by calculating the energy of an idealized membrane hole. The relevant dimensionless 
parameter P is the ratio of the line tension to the electrostatic energy per screening length of edge P — e w n 2 ^/(j 2 . We 
show that when line tension dominates, holes close. When P< 1, large holes grow. At intermediate P w 2 metastable 
holes exist due to the competition between energies. We then discuss the stability of a straight membrane edge to 
sinusoidal perturbations with wavenumber q. We show that at small P, and in the long wavelength limit q <C n, small 
perturbations lower the overall energy. In this regime, a flat membrane with a free edge will generically break into 
pieces of minimum size 0(k _1 ). The final section of the paper relates these results to the ghost experiments. 

To compute the electrostatic energy of a hole, we use the linearized Poisson-Boltzmann equation jO] V 2 </> = k 2 4>. 
The boundary conditions are <f> — > as z — > 00, d z <j)(z = 0) = —2ir<j/e w , r > R and d z (f>(z = 0) = for r < R. The 
latter comes from symmetry of the potential about 2 = 0. We Hankel transform to find 



where Jo and J\ are Bessel functions. Without a hole (R = 0), the potential reduces to that of an infinite charged 
sheet. 
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The electrostatic energy is U — 2tt f °° rdr a(f>(z = 0); we assume the total charge and the membrane area remain 
constant. The electrostatic contribution competes with the energy from line tension, 2it"fR. The energy change, 
nondimensionalized by 2ira 2 /(e w K 3 ), depends on the parameter P = e w K 2 j/a 2 and the dimensionless radius £ = kR 

AU = P£- + 2k? r dX f\ {x) ] 2 (2) 
Jo x^x z + 1, 

When the hole is small (£ <C 1), the integral in the third term is approximately 4/(37r), so the energy change of a hole 
is AU « P£ — 7r£ 2 + 8/3£ 3 . Thus, for small holes the line tension dominates over electrostatic effects. For larger £, the 
energy is shown in Fig. 3. When P > 2.2, the energy increases monotonically as a function of £, so holes always close. 
For 2.0 < P < 2.2, we find a local energy minimum with £ ~ 1. In principle, holes in a membrane could be observed 
in this region of metastability. In practice, however, we believe that the range of parameters in which such holes exist 
is prohibitively narrow. When P < 2, the energy decreases linearly at large £, so that large holes will grow. In this 
regime there is an activation energy barrier to hole formation of size e w K^ 2 a^ 2 /2. 
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FIG. 3. Energy of a hole as a function of dimensionless radius k,R, for P = 1.8, 2, 2.05 and 2.25. A typical barrier height is 
10 kT (for P — 1, which corresponds to n=l mM, 7 = 10 -6 erg, and a = 0.2 e/nm 2 ). 

A membrane with a growing hole ultimately has a long edge; is this edge stable? We consider a straight edge 
segment (at x = 0) modulated by a sinusoidal perturbation x — ecos(qy). For this problem the charge density 
is a — <Jo<d(x + ecosgy), where is a step function. The energy follows from the linearized Poisson-Boltzmann 
approximation by expanding the potential and the charge density in powers of e. The first nonzero correction to the 
energy from the edge perturbation occurs at 0(e 2 ), and this change in the electrostatic energy (in dimensionless units) 
due to the edge modulation is — e 2 K 3 g _1 ln(l + n~ 2 q 2 ). 

The perturbation increases the length of edge by AL = ir/2e 2 q, resulting in an energy increase 7AL from line 
tension. For large q line tension dominates, increasing the energy. For small q -C k, the total energy change is 

AU = ^f(^-?j + 0(q 2 ), (3) 

If (a) electrostatic effects dominate (P < 2) and (b) the perturbation wavelength is greater than the screening length, 
modulations of the edge grow. The growth continues as long as the circumference of membrane fragments is longer 
than the screening length. When the fragment radius is of order 27tk~ 1 the instability will stop, and the fragment may 
close to form a vesicle. In the absence of charge, the fragment size must be > 10 nm jlj|, so the energy associated 
with the free edge is larger than the bending energy required to close into a vesicle. With electrostatic effects included, 
it seems additional forces are necessary to promote vesicle formation. A candidate for red blood cells is spontaneous 
curvature: the higher charge density on the inner face of the bilayer || means the membrane prefers to bend towards 
the outer face B. Understanding the final vesicle sizes requires analyzing both the electrostatic edge instability 
and spontaneous curvature — in the Lew experiments, || spontaneous curvature effects play an important role in the 
experimental phenomenology. 
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Observations of stable holes in red blood cell ghosts motivated this work |lfij]. Human red cells, if burst by osmotic 
stress and placed in a low-salt buffer, break into vesicles 0, which typically show the cytoplasmic (inner) side of the 
membrane facing outwards. Studies of this process illustrate that vesiculation is preceded by the formation of 
stable holes in the membrane, whose size depends on the salt concentration of the solution and the membrane charge 
density. Lieber and Steck |8| showed that (a) the hole size increases with decreasing salt concentration; and (b) a 
decrease (increase) in the effective membrane charge density — via charged intercalators inserted into the membrane — 
decreases (increases) the hole size. Our model qualitatively explains these trends, but uncertainties in our model 
and in experimentally measured parameters prevent precise numerical calculation. We have neglected (a) differences 
in charge density on the two sides of the membrane, (b) the insulating interior of the bilayer and (c) nonlinear 
effects. Although there are experimental indications that some of these effects might be relevant jT^], in the absence 
of quantitative information for the line tension and charge density, the simple model presented here is adequate. 

Fig. 4 shows the stability boundaries as a function of the surface charge density a and the salt concentration n. We 
assume 7 = 10~ 6 dyn, consistent with recent experiments [ p"3| |. The phase boundaries follow a oc n 1 / 2 . For typical 
membrane charge density a ~ 0.2 e/nm 2 , the crossover point P ~ 2 occurs at a salt concentration of 2 mM, which is 
in the range of the Lieber-Steck experiments. When the energy barrier for hole formation is of order a thermal energy 
ksT, thermal fluctuations produce holes in the membrane. The dashed line in Fig. 4 shows the borderline for this 
instability |Q . Since the barrier scales like 7 2 , the exact position of this important borderline is sensitive to the line 
tension 7. 
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FIG. 4. Phase Diagram, depicting regimes of shrinking, growing and metastable holes, with 7 = 10 dyn. The upper 
and lower solid lines depict P — 2 and 2.2. The electrostatic edge instability is present everywhere above the solid lines. 
Note that the salt concentration (assuming a 1:1 salt such as NaCl) is related to the Debye length via k 2 = n/£ 2 , where 
I 2 = e w k B T/(8ne 2 N), with N Avogadro's number. 

The electrostatic model qualitatively explains the dependence on salt concentration and membrane charge density 
in the experiments. However, electrostatic effects arc not sufficient: the edge instability predicts that holes either 
grow indefinitely (resulting in vesiculation) or close completely. The observations of stable holes for long time 
periods contradict this prediction and indicate that a stabilizing effect is necessary. 

Several aspects of the experiments indicate that the spectrin skeleton — a protein mesh anchored to the membrane — 
is this stabilizing element. The electron micrographs of Lew et. al. |)| show that the spectrin network is intact before 
vesiculation occurs, and as the ghosts disintegrate into vesicles the spectrin mesh detaches from the membrane. 
Lieber and Steck || found that proteins known to covalently cross-link spectrin stabilize the hole, preventing changes 
in radius as the ionic strength of the solution is changed. Conditions known to promote breakdown of spectrin (such 
as increase in temperature or digestion by enzymes Q) cause the ghosts to vesiculate. If the spectrin provides the 
restoring force needed to stabilize the membrane, force balance implies k s R = —dll/dr, where k s is the spectrin area 
expansion modulus and U the hole energy computed above. At small P, we find 

47TCT 2 ^ 2 1 

R =^1J_L. (4) 
Figure|| shows Lieber and Steck's data ||], replotted on a double logarithmic scale. 
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FIG. 5. Radius of hole versus salt concentration from Ref. 0. The solid line is R ~ n _1 . 

The data obey this law over a decade in hole size. To check the consistency we fit the prefactor of the scaling law. 
A best fit gives the prefactor c = 7.0 • 10~ 6 cm M. For a in units of e/nm 2 and k s in dyn, we find a 1 jk s — 0.02. This 
ratio is consistent with the charge density a = 10~ 2 and bulk modulus k s = 10~ 2 from independent measurements 




The mechanism of vesiculation by electrostatic edge instability does not require spontaneous curvature; hence it 
differs from other explanations of vesiculation There are conditions under which vesiculation appears to be 

driven by spontaneous curvature JTsj ] . It seems likely that both mechanisms for vesiculation operate in practice. 

For physiological conditions, membrane charge densities range from —0.03 to —0.24 |e|/nm yjMpl and NaCl 
concentrations are around 150mM. The phase diagram Fig. 4 predicts that such membranes should be stable, although 
conditions in vivo are within an order of magnitude of the stability boundary. This indicates that the edge instability 
could be important in living systems. Partition of the nuclear envelope (the membrane surrounding the nucleus) during 
cell division provides a possible example [ p2|j23| . Because it is a double bilayer IJ, its effective line tension is much 
smaller than a typical membrane. The qualitative behavior of the nuclear envelope is similar to the phenomenology 
described here: during most of the cell cycle, the nuclear envelope is attached to a protein meshwork called the nuclear 
lamina |2^| . Near the start of cell division, the lamina disassembles — and the nuclear envelope is observed to fill with 
holes and vesiculate p3[ . The many vesicles are divided between the two daughter cells during division. 
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